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Abstract 



We discuss the influence of Dirac sheets and zero-momentum modes on the 
gauge variant photon correlators Y{r;p) with p ^ and p = in the pure gauge 
{7(1) theory. A special attention has been paid to the (3- and volume-dependence 
of this influence. Numerical simulations are performed on 12 x 6 3 and 24 x 12 3 
lattices at different /3's in the Coulomb phase. 

b . 

1 Introduction 

Lattice approach |] gives a possibility to calculate gauge invariant objects without 
gauge fixing. However, numerical calculations of gauge dependent objects, e.g., 
gauge field correlators T(r; p) , can also be of interest. For example, the comparison 
of gauge variant objects with that calculated perturbatively can provide some deeper 
insight into the structure of the lattice theory. In particular, the role of lattice 
artifacts and/or Gribov copies |2|] can be investigated. 

Compact i7(l) pure gauge theory in the Coulomb phase provides a unique 'test 
ground' for the lattice approach. Indeed, in the weak coupling limit one expects to 
find a trivial theory of the free non-interacting photons. Therefore, at sufficiently 
small coupling a gauge variant transverse correlator T T (r]p) is supposed to fit the 



perturbative expression given in eq. ( |3.5| ). 

A few years ago it has been shown |J that in the Coulomb phase some of the 
gauge copies produce a photon correlator IY(r;p ^ 0) with a decay behavior 



1 



inconsistent with perturbation theory. Numerical study |4| has shown that there 
is a connection between 'bad' gauge copies and the appearance of configurations 
with periodically closed double Dirac sheets (DDS). The explanation of this effect 
0] is connected with a nontrivial 'vacuum' structure of the classical compact U(l) 
theory, i.e. with the existence of the nontrivial stable classical solutions ('vacua'), 
some of these solutions - DDS - being Gribov copies of the trivial 'vacuum' with 
zero field. 

Another interesting observation is connected with the r-dependence of the zero- 
momentum gauge variant correlators T(r) . The analysis of the r-dependence of 
the correlator T(r) in the SU(3) lattice gauge theory drove the authors of || to 
the conclusion that T(r) is consistent with the propagation of a massive gluon. On 
the other side, in the U(l) theory the correlator T(r) exhibits a similar behaviour 
0. However, one can hardly expect the appearance of the massive photon in the 
Coulomb phase in the U(l) theory. It has been shown in ref. |5| that (at least, 
in the case of the U(l) theory) the r-dependence of the correlator T(r) can 
be explained taking into account zero-momentum mode of the gauge field. The 
influence of the zero-momentum modes on the gluon propagators in nonabelian 
theories has been pointed in || [LC], [TT[ . 



In this paper we discuss the effects connected with Dirac sheets and zero- 
momentum modes in more details. The main questions addressed in this paper 
are the following. 

- What is the /^-dependence of the frequency of the Dirac sheets appearance ? 
Does the contribution of the Dirac sheets to the gauge variant correlator T(r;p) 
become more (or less) important when the coupling (3 is increased? 

- What is the volume dependence of the correlator IY(r;p) ? Does the influence 
of DDS become smaller in the thermodynamical limit ? 

- What is the influence of zero-momentum modes on the zero-momentum correlator 
T(r) ? How does this influence depend on (3 and lattice size ? 

Throughout this paper a 4d lattice with periodic boundary conditions is con- 
sidered. Np is the lattice size in the direction /i, and V4 = N1N2N3N4. The lattice 
derivatives are d fi f(x) = f(x + jl) — f(x) and d^f(x) = f(x) — f(x — ft), and the 
lattice spacing is chosen to be unity. 

2 The action and classical solutions 

The standard Wilson action S(U) is 

5(f/) = ^EE( 1 - cos M , (2- 1 ) 

X fJ,>U 
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where U Xfl = e l6xti G Z7(l) are link variables, O x .^ v = d^9 xu — d u 9 x ^ are the plaquette 
angles and j3 = 1/g 2 ■ This action is the part of the full QED action Sqed, which 
is supposed to be compact if we consider QED as arising from a subgroup of a 
non-abelian (e.g., grand unified) gauge theory flTjJ . 

The plaquette angle P = x . ifiu can be split up: P = [Op] + 2irnp, where 
[6p\ G (— 7r; 7r] and rip = 0, ±1, ±2. The plaquettes with rip ^ are called Dirac 
plaquettes. The dual integer valued plaquettes m x ^ u = \e^ U p a n x pa form Dirac 
sheets @. 

In a perturbative language the existence of the gauge copies can be interpreted 
as a problem of the gauge copies of the 'vacuum' configurations, i.e. the solutions of 
the classical equations of motion. Assuming that every configuration is some small 
fluctuation about the corresponding 'vacuum' one can find all gauge copies of this 
configuration provided all gauge copies of the 'vacuum' are known. 

The classical equations of motion are 

$>sin^ = . (2.2) 

V 

Evidently, zero-momentum modes = M are the solutions of these equations. 

Choosing the Lorentz (or Landau) gauge YLp d^0 Xil = , one can write expicitely 
a solution corresponding to a single Dirac sheet || 



where /C M = 2 sin ^ and K? ± = K\ + K 2 . The two-dimensional vector R = (Ri, R 2 ) 
corresponds to the position of the Dirac plaquette in the (x±] x 2 ) plane: 

^12 = 2tt • - A ; t " 



N X N 2 



Of course, 0' = 4>u + 0^ is also a solution of eq's. 



The single Dirac sheet solution 0^ (R) corresponds to a local minimum of the 
action, i.e. that it is stable with respect to small fluctuations. The existence of the 
long-living metastable states corresponding to single Dirac sheets was observed in 
simulations in the pure gauge U(l) theory |T4|. 



The classical gauge action is S(0 d ) = fl — cos A^ . On a symmetric lattice 
Ni — N 2 — N 3 — — > oo the action S(0 cl ) is non-zero and finite : 
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4tt 2 

S(9 d ) = — < oo . 
T 

Double Dirac sheet solution of the classical equation of motion consists of the 
two single Dirac sheets with an opposite orientation of the flux : 

9*(R a ; R») = e:[(R a ) - 9*(R b ), i = l-2, (2.5) 

where the vectors R a and R b correspond to the two Dirac plaquettes in the plane 
{x\; x 2 ). It is easy to see that 

9* 12 (R a ; R b ) = 2tt ■ [8 3x . &a - 5 3x . Rb ] . (2.6) 

The double Dirac sheet 9^. i (R a ;R b ) has a zero action: S(6 cl ) = 0. 

Gauge transformations can shift the Dirac sheets and change their form. For 
example, the 'big' gauge transformation function Q x 

Q x = — y (l-e- iq A ; (2.7) 

shifts the single Dirac sheet in the Xi-direction. The gauge transformation Q x in 
eq. (|2.7| ) applied to the zero-field creates a double Dirac sheet as in eq. fl2.5j ) with 
R« = R a ndR b = R-i. Therefore, fl* (#; R b ) is a Gribov copy of the zero 
solution 9^ = 0. It is not difficult now to obtain general Dirac sheet solutions, i.e. 
the Dirac sheets curved in the four-dimensional space ||. 

3 Tt(t\p) and Dirac sheets 

In lattice calculations the usual choice of the Lorentz (or Landau) gauge is 

4 

53 d» sin 9 xfM = 0, (3.1) 
which is equivalent to finding extremum of the functional F(9) 

F ( d ) = TfT, F M ; F M = ~E[ C0S ^ + C0S ^-aJ (3.2) 

with respect to gauge transformations U X}X — > A x U x ^A x+fl ; A x = exp{iQ x } . It 
is important to note that all link angles 9 Xfl are compact variables (— 7r < 6 XfJi < tt), 
and the gauge transformations 8 Xfl — > 9 X[X — dyVL x are understood modulo 2ir. 

To find extremum of F{9) we carried out a standard overrelaxation gauge 
cooling procedure with parameter a tuned to minimize the number of the cooling 
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steps for given (5 and volume. In our computations values of a were chosen to 
be between a = 1.66 and a = 1.84 . The stopping criterion was 



max 



dfi sin 9 : 



Xfl 



< 10~ 5 and — Y 



The photon correlator r^(r;p) is 



^ sin 9. 



xp 



< 10" 



where t © r = (t + r) mod A/4 and 



AZ 



1 N4-1 

Eto*©r;p)0 M (t;p)); (3.3) 



4 t=0 



e^(r;p) = ^ e -^-3Pf . s in# 



XfJ, ; 



/i = 1,2,3 



(3.4) 



and, clearly, (O m ) = 0. 

Let us choose the momentum p = (0;p 2 ;0) with p 2 7^ 0, and // = 1. The 
perturbative expansion about the zero solution of the classical equation of motion 
®xu = O5 i- e - the standard perturbation theory, gives in the lowest approximation 



(3.5) 



(3.6) 



rf T \r- p) ~ e" r ^ + e" (iV4 - T) ^ , 
where the energy E p satisfies the lattice dispersion relation 

smn — = > sin — . 

2^2 

z i=i z 

It is easy to see that the expansion about a Dirac sheet solution gives a contri- 
bution to the correlator very different from that in eq. ( |3.5| ). 

As an example, let us choose the flat double Dirac sheet with the space-like 
Dirac plaquettes 9^(Ri, R 2 ) as defined in eq. fl2.5|). In this case the correlator in 
the double Dirac sheet background is 



3 g± 



(3.7) 



where 



m = ivfll ? sin * ; *® = a^ ? e " rx ^ cos * ' (3 ' 8) 

and 
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G(t; q) 
Go(r;q) 



1+ /C?(g) d 



2sinhi? g dE q ^ 



2sinh£ q 1 - e- N * E * 



G (r;q); (3.9) 



_|_ e -(N4-r)E g 



+ ... . 



This expression is obviously different from that given in eq. ( p.5|) . In particular 



the first term in the r.h.s. in eq.(3.7) produces the 'shift upward' of the correlator 
with respect to the perturbative expression. 

In our computations we monitored the total number of the Dirac plaquettes 
N$$ for every plane (fiu) and Np>p = max( p ) N^p . By a detailed investi- 
gation of the number and the location of Dirac plaquettes in the Coulomb phase, 
we observed pairs of Dirac sheets closed by periodic boundary conditions. In the 
extreme case of a minimal surface this is a pair of periodic Dirac sheets occupying 
parallel planes and having corresponding Dirac plaquettes with opposite signs. The 
appearance of periodically closed double Dirac sheet means that the number of the 
Dirac plaquettes in a given plane (fiu) has to be 



Nnp > 



DP ~ N,N U 

For example, on a 12 x 6 3 lattice, the minimal number of Dirac plaquettes is 
JVj™ n = 72 if the Dirac sheets are timelike ('short' DDS), and iVj™ n = 144 if the 
Dirac sheets are spacelike ('long' DDS). 

Figure 0a shows a time history of the number of the Dirac plaquettes N^p on 
the 12 x 6 3 lattice at (3 = 1.1 . For every configuration the number Np,p has 
been measured after the gauge fixing procedure. The number of configurations with 
Dirac sheets (N DP > 72) is ~ 20% . Long DDS (N DP > 144) are not as frequent 
as short DDS. However, as we shall see later the influence of the long DDS on the 
photon correlator is much stronger. 

To convince ourselves that the appearance of the Dirac sheets is not an artifact 
of the small volume and/or small (3- values we repeated these calculations on the 
lattices 12 x 6 3 and 24 x 12 3 at /3-values between (3 = 1.1 and (3 = 10 . We 
observed that the probability to find DDS shows rather weak dependence on the 
coupling. We observed also that the probability to find DDS increases slightly with 
increasing of the volume (e.g., at (3 = 2 the number of configurations with DDS 
is ~ 17% on the 12 x 6 3 lattice and ~ 25% on the 24 x 12 3 lattice). In 
Figure |I|b we show a time history of the number of the Dirac plaquettes Np,p on 
the 24 x 12 3 lattice at (3 = 2. This time history looks less 'noisy' as compared 
with that in Figure ffk because there are less monopoles at (3 = 2. 
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In Figure |2]a we show the transverse correlator T(r;p) with p = (0; |^;0) 
at [3 — 1.1 on the 12 x 6 3 lattice. Solid line corresponds to the lowest order 
perturbative expression given in eq.( |3.5| ). Averaging over all configurations (filled 
circles) gives noticeable disagreement with the free photon case. On the other hand, 
excluding from the statistics configurations containing Dirac sheets gives an aver- 
age correlator (opened circles) which fits nicely the zero-mass photon correlator. 
Squares represent a correlator obtained by averaging over configurations contain- 
ing long Dirac sheets. The destructive influence of the Dirac sheets exhibits little 
changes with increasing (3 and the lattice size. As an example we show in Figure 
|2|b the transverse correlator T(r; p) at (3 = 2 on the 24 x 12 3 lattice. 

The results presented in this section could be summarised by saying that the 
contribution of configurations with DDS 'spoils' the photon correlator and leads to 
a wrong dispersion relation inconsistent with the dispersion relation for the massless 
photon. Dirac sheets represent an example of lattice artifacts which can lead to a 
wrong interpretation of the results of numerical calculations. 

4 Zero— momentum correlator 

The (connected) zero-momentum correlator T(r) is 

r„(r) = (0^)0,(0)) ; 0,{t) =J2sin6 x , , (4.1) 

X 

where r is the separation in one of the four euclidian directions and x corresponds 
to the three complementary directions. 

A usual way to define effective masses m e //(r) is 



cosh 


m e// (r)(r + l-|7V 4 )" 


r„(r + i) 


cosh 


m eff(r)(T- iAT 4 )_ 





In the SU(3) theory the large-r behaviour of m e //(r) is supposed to describe the 
'gluon mass' m g if r is chosen along the 'temperature' direction x 4 and \i = 1, 2, 3 

§• 

The definition of the effective mass in eq. (^.2|) is such that for a correlator 
behaving like T^r) ~ e~ Mr + e ~ M ( N ^- T ) ( one obtains m e ^ = M . However, on the 
periodic torus the gauge variant correlator Tj(r) cannot be represented as a sum 
of exponents. Zero-momentum modes make their contribution to this correlator 
adding a nonzero positive constant as in eq. (|4.7|) 0. As we shall see later this 
contribution cannot be discarded. 

First, it is rather instructive to look at the distributions of the operators 0^(t) 
defined in eq.( |4.1| ) and the correlations between different time-slices r, where the 
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index k stands for the k th measurement. The width of the distribution of the zero- 
momentum operators is not small, and shows very weak dependence on the coupling 
constant in the interval between (3 = 1.1 and (3 = 10.0. As an illustration we show 
in Figure |3] a 'scatter-plot' of the operators 0^(t) with r = 1 and r = 2 at 
(3 = 2 on a 24 x 12 3 lattice. For every measurement /c the operators QV*\t) can 
be represented as 

{k \r) = C {k) +50 {k) (r) , (4.3) 

where 50^ k \r) are small fluctuations about the constant contribution . The 
width of the stripe in Figure |3| reflects the amplitude of these fluctuations. 

It is worthwhile to note that correlations between operators 0^ k \r;p) with 
nonzero momentum are very different from that for 0"(t) , namely the scatter 
plot for 0«(r;p) consists of a small 'blob' at the center of coordinates (is not 
shown) . 

It is not difficult to calculate the gauge variant correlator in the lowest approx- 
imation. The average of any functional $(#) is defined as 

($) = ^ f[d6 x »] $(9) ■ e- Seffie) , (4-4) 

where S eff = S + S gf + S FP , S gf = -±p E*(£/j % sin^) 2 , and S FP is the 
contribution from the corresponding Faddeev-Popov determinant. The choice a = 
corresponds to the Lorentz gauge. 

Lorentz gauge condition, does not exclude the appearance of the zero-momentum 
modes in the periodic volume. To calculate ($) perturbatively one should keep the 
zero-momentum mode under control in the perturbation expansion. This can be 
achieved by repeating the Faddeev-Popov trick 
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r* 4 f 1 

/ n<% ex p^ — n 

J -%=1 l 6 A 4 



(4.5) 



and making the change of variables 9 XI1 = M + gA Xfl . The average functional ($} 
can be now represented in the form 

<*> ~ f [<%] f[dA x ,} Y[S(Y / A x M ( p + gA)e~ s ^ + ^ . (4.6) 

In eq. ( |4.6|) one can safely expand in powers of g 2 and extend the limits of integration 
of A Xfl to [— oo; oo]. The zero-momentum modes are not gaussian (the gauge action 
S(9) does not depend on M ), and the integration over M in eq.( [4.6|) should stay 
compact. 

Applying a perturbative expansion about the constant (zero-momentum) modes 
M to calculate the zero-momentum correlators r^(r) one obtains 



S 



r M (r) ~ (sin 2 M > ■ (1 - b,g 2 ) + |^(cos 2 M > • T'^r) , (4.7) 



where 



T'(t) 



9 2 Vi 



^3 

2V A 



^ 4sin 2 2f 



z = l;2;3; 







P4^0 

at a 



0. 



(4- 



and V3 = N1N2N3. Coefficients represent tadpole contributions : 



2V 4 



E 



1 



JC,, = 2 sin 



(4.9) 



The coefficients (sin 2 0^) and (cos 2 (j)^) have been numerically calculated from 
distributions of the zero-momentum operators. 

As an example we show in Figure |] the correlator T(r) = Tj(r) at (3 = 2 on 
a 12 x 6 3 lattice. Averaging over configurations without DDS (filled circles) gives 
a nice agreement with correlator defined in eq.( |4.7| ) (solid line). The agreement 
becomes even better with increasing (3 . The increasing of the lattice size does 
change this conclusion. It is worthwhile to stress that this agreement is impossible 
without taking into account zero-momentum modes. Correlators defined on the 
configurations with DDS (diamonds and triangles) are very different from that 
defined on the configurations with DDS excluded. 



5 Conclusions 

We have studied the influence of Dirac sheets and zero-momentum modes on the 
gauge variant correlators T(t;p) with zero and nonzero momenta. All calculations 
have been performed on 12 x 6 3 and 24 x 12 3 lattices at different (3 's in the 
Coulomb phase. 

We confirm our previous conclusion ||, [J about the role of the Dirac sheets : 
gauge copies which possess double Dirac sheets contribute to a wrong ('nonphysical') 
behaviour of the gauge variant photon propagators T(r;p) at p = and p 7^ . 

We have checked that with the increasing of the coupling from (1 — 1.1 up to 
j3 — 10 on a 12 x 6 3 lattice and from /3 = 1.1 up to j3 = 2 on a 24 x 12 3 
lattice the frequency of DDS does not change significantly. This frequency grows 
somewhat with the increasing of the volume. Because of configurations with DDS 
the numerically measured gauge variant correlators do not correspond to the free 
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massless photons. Therefore, to avoid nonphysical conclusions one should exclude 
configurations with Dirac sheets out of consideration. 

The contribution of the zero-momentum modes is of crucial importance for the 
interpretation of the zero-momentum correlators T A1 (r) . These modes cannot 
be discarded at all (3's and volumes we employed. A formal application of the 
definition of the effective mass m e // given in eq.( |Q|) can produce a misleading 
interpretation. 
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Figure captions 

Figure |l|. Time history of N DP at f3 — 1.1 on the 12 x 6 3 lattice (a) and at (3 — 2 
on the 24 x 12 3 lattice (b). 

Figure ||. T t (t;p) at (3 = 1.1 on the 12 x 6 3 lattice (a) and at (3 = 2 on the 
24 x 12 3 lattice (b). 

Figure [| Scatter-plot for operators 0{t = 1) and 0{j = 2) at (3 = 2 on the 
24 x 12 3 lattice. 

Figure |. T(t) at (3 = 2 on the 12 x 6 3 
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Figure 1: Time history of N DP at (3 — 1.1 on the 12 x 6 3 lattice (a) and at (3 — 2 
on the 24 x 12 3 lattice (b). 
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Figure 2: T t (t;p) at (3 = 1.1 on the 12 x 6 3 lattice (a) and at /3 = 2 on the 24 x 12 3 
lattice (b). 
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Figure 3: Scatter-plot for operators 0(t = 1) and 0(t = 2) at (3 = 2 on the 
24 x 12 3 lattice. 
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Figure 4: T(r) at (3 = 2 on the 12 x 6 3 lattice. 
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